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Abstract 

We explore the consequences of multi-trace deformations in applications of gauge- 
gravity duality to condensed matter physics. We find that they introduce a powerful 
new "knob" that can implement spontaneous symmetry breaking, and can be used to 
construct a new type of holographic superconductor. This knob can be tuned to drive 
the critical temperature to zero, leading to a new quantum critical point. We calculate 
nontrivial critical exponents, and show that fluctuations of the order parameter are 
'locally' quantum critical in the disordered phase. Most notably the dynamical critical 
exponent is determined by the dimension of an operator at the critical point. We 
argue that the results are robust against quantum corrections and discuss various 
generalizations. 
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1 Introduction 

Over the past couple of years, gauge/gravity duality has been applied to a number of prob- 
lems in condensed matter physics (for reviews see PQ |2J |3]). An important feature of some 
condensed matter systems is the existence of quantum critical points, marking continuous 
phase transitions at zero temperature. One goal of the present work is to introduce and study 
a new mechanism for generating quantum critical points in the context of gauge/gravity du- 
ality. We will see that the behavior near the critical points is described by nontrivial critical 
exponents and goes beyond the usual Landau-Ginzburg- Wilson description of phase transi- 
tions at zero temperature. 

A second goal is to introduce a new type of holographic superconductor. The key ingre- 
dient in constructing a gravitational dual of a superconductor is to find an instability which 
breaks a U(l) symmetry at low temperature and causes a condensate to form. Previous 
constructions have started with a charged anti de Sitter (AdS) black hole which has such an 
instability when coupled, e.g., to a charged scalar field jH El E]. We will show that there is 
another source of instability which applies even for Schwarzschild AdS black holes. So these 
new superconductors can exist even at zero chemical potential and no net charge density. 

Both of these goals are achieved by adding a multi-trace operator to the dual field theory 
action^] For example, given a (single trace) scalar operator O of dimension A_ < 3/2 in a 

2 + 1 dimensional field theory (which will be our main focus) one can modify the action 



where for convenience later the coupling will be rescaled as R = 2(3 — 2A_)/t. Since this 
is a relevant deformation, it is unnatural to exclude such a term, and it has important 
consequences. If O is the operator dual to the bulk charged scalar field in conventional 
holographic superconductors, then adding this term (with k > 0) makes it harder to form 
the condensate and lowers the critical temperature. We will see that in some cases T c vanishes 
at a finite value of k = k c . This defines a new quantum critical point which we will study in 

1 For another recent discussion of multi-trace operators in gauge/gravity duality, see [7J. 
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detail. Since T c can be quite large at K — 0, adding this double trace perturbation introduces 
a sensitive new knob for adjusting the critical temperature. 

For k > k c (and nonzero chemical potential //), the ground state is described by the 
extremal Reissner-Nordstrom (RN) AdS black hole which has an emergent AdS2 geometry 
in the IR. For k < k c there are various possible IR geometries depending on details of the 
bulk potential. However, as k approaches n c from below, the bulk solution develops an 
intermediate AdS2 geometry. It is this intermediate region which controls the behavior near 
the critical point. For example, we will show that the critical exponents do not take mean 
field values, but are determined by the scaling dimension of certain operators in the + 1 
dimensional CFT dual to this region. This is closely analogous to the way properties of 
holographic non- Fermi liquids [8j [9j HO] were described in terms of a dual + 1 dimensional 
CFT [TTJ [12]. In addition, the instability for k < k c can be interpreted as turning on a 
double trace term with negative coefficient in the + 1 dimensional CFT dual to this region. 

Let us contrast this with the usual argument for why the RN AdS black hole becomes 
unstable at low temperature in the presence of a scalar field [13J. In the AdS2 near horizon 
geometry of the T = solution, the scalar field has an effective mass m eS which depends 
on the original mass m and charge q of the scalar field. When this effective mass squared 
is below the Breitenlohner-Freedman (BF) bound [H] for AdS2, this near horizon region 
becomes unstable. Since m 2 is above the BF bound of the asymptotic AdS^ geometry, the 
asymptotic region is stable, and the solution settles down to a black hole with scalar hair. 
It is now clear that this argument is sufficient but not necessary. It overlooks the possibility 
of instabilities with above the BF bound which are allowed due to modified boundary 
conditions for the scalar field. The boundary conditions may be modified due to the addition 
of a multi-trace deformation in the dual field theory [T5l IT6] , or simply due to alternative 
quantization of the bulk theory [PT] . 

It was widely believed that if one added a double trace term with k < 0, then the 
theory would not have a stable ground state. However, we have recently shown that this 
is not necessarily the case [IB]. For a large class of dual gravity theories, there is still a 
stable ground state with (O) ^ when the boundary conditions correspond to n < 0. 
As one increases the temperature, there is a second order phase transition to a state with 
(O) = 0. This provides a new mechanism for spontaneously breaking a U(l) symmetry and 
constructing novel holographic superconductors. This mechanism does not require a charged 
black hole and works for Schwarzschild AdS as well. In other words, one can set /i = and 
still break the U(l) symmetry at low temperature. The critical temperature is now set by 
Kj. We will discuss some properties of these novel holographic superconductors in section 2. 
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It is worth pointing out that the coupling k, as the coefficient of the square of the order 
parameter, is the usual tuning parameter in the context of Landau-Ginzburg theory. Also 
if O is a gauge invariant trace of a fermion bilinear then the double trace is a 4 fermion 
interaction, a natural interaction to consider. 

Although our discussion so far has focussed on the case where the operator O is charged, 
our results apply equally well when O is neutral. In this case, the ordered phase breaks 
a Z 2 symmetry. More generally, one can imagine different symmetry breaking scenarios 
where for example O could be part of a triplet of operators forming a representation of 
SU(2) which is spontaneously broken at low temperature. This is a particular attractive 
possibility as outlined in [191 120], since in many condensed matter systems 577(2) spin is a 
global symmetry (ignoring spin orbit effects.) Including an exact global SU(2) symmetry 
then allows us to model magnetism in a holographic setup. The triplet in which one embeds 
O can be interpreted as the staggered order parameter associated with anti-ferromagnetic 
transitions. Since the boundary theory has a global SU(2) symmetry the bulk will have an 
SU(2) gauge symmetry distinct from the £7(1) electromagnetic charge. For the rest of this 
paper the SU{2) gauge fields and triplet structure of the order parameter will not play a 
roll. 

It is useful to bear this possibility in mind, in particular so we can compare our results 
to quantum phase transitions in metallic systems, where anti-ferromagnetic order plays an 
important role. The "standard" theory of which was given in [211 E2] is based on the renor- 
malization group Landau-Ginzburg paradigm. However experimental measurements (see for 
example [231 EI], and references therein) of heavy fermion systems with quantum critical 
points show that the "standard" theory can break down, as a result new theoretical meth- 
ods are required. Subsequently several different methods were developed (see for example 
[251 12"6l [271 128].) One such method [251 EE] which is formally justified in a large d expansion 
[29] shows local quantum critical behavior similar to the new quantum critical point that we 
find. It will be useful to compare and contrast our results to those of the standard theory 
and the other theoretical methods used in the study of quantum criticality in heavy fermion 
systems. 

Since the finite density normal phase we consider is governed by AdS 2 x R 2 in the IR, 
at zero temperature the theory has a finite entropy density. This has lead many people 
to suggest that this state must not be the true ground state, since otherwise one finds 
unnatural violations of the third law of thermodynamics. Of course this may be natural in 
the context of applications to heavy fermion system where superconductivity instabilities can 
be observed close to criticality. So it may be that the state we work in is the correct one for 
a large range of temperatures, but ultimately at low temperatures something else takes over. 
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Indeed as is discussed and extended in this paper AdS2 has many possible forms of instability. 
This however motivates us to attempt to extend our results in various directions to directly 
address this problem. One extension we consider is adding a magnetic field. We show that 
while a magnetic field suppresses the superconducting instability, it can enhance the neutral 
(anti-ferromagnetic) instability. Another extension we consider is replacing AdS2 x R 2 with 
other possible IR geometries, such as a Lifshitz geometry which does not have finite entropy 
density at zero temperature. We find our results are rather robust here. 

The organization of the paper is as follows: In the next section we show how double trace 
deformations can induce spontaneous symmetry breaking and use this to construct a new 
type of holographic superconductor with zero net charge density. We then extend this to the 
finite density case and show that the coefficient of the double trace deformation provides a 
sensitive knob by which one can tune the critical temperature T c to zero. In section 3 we 
study the new quantum critical point that arises and analytically compute the nontrivial 
critical exponents. In section 4 we numerically construct the backreacted geometries that 
correspond to the ordered (condensed) phase away from the phase transition. We confirm the 
critical exponents near the critical point. In the discussion section, we summarize our results 
and discuss generalizations to magnetic fields and Lifshitz normal phases. The Appendices 
contain additional details. 

2 Double trace deformations 

In this section we begin by emphasizing the simple under appreciated fact that double trace 
deformations are useful for studying symmetry breaking in gauge/gravity duality. We then 
note that this system provides a simple holographic model for superconductivity with zero 
total charge density. With nonzero charge density, we show that double trace deformations 
introduce a new parameter by which one can tune the critical temperature of the supercon- 
ductor. 

2.1 Gravity setup and boundary conditions 

The theory we study is gravity in 3 + 1 dimensions with a negative cosmological constant, 
a U(l) gauge field, and a scalar field \l/ which may or may not be charged under the U(l) 
symmetry. By general arguments of gauge/gravity duality this theory is dual to a CFT 2+ \ 
with a conserved current operator J M and a scalar operator O. 

The action is that of the Einstein-Abelian Higgs model with a negative cosmological 
constant, where we parameterize the phase and modulus of the charged scalar as ^> = ipe %e , 
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following e.g. [30] 



S = J d'x^g (r - \g^)F 2 - (V^) 2 - J^)(Ve - qAf - V(rl>)\ (2.1) 

We require that the coupling functions G and J and the potential V be even functions of ip, 
since we need to preserve our U(l) symmetry. We will assume an expansion of the form 

l/ = -6 + mV + 0(^), G= l + # 2 + 0(V> 4 ), J = ^ 2 + C(^ 4 ) (2.2) 

where we have set the AdS± radius to one. The coefficient of if) 2 in J is fixed by regularity 
at ijj = 0. This is all that we will need to determine the behavior near the critical point in 
section 3. We will specify the potential and coupling functions more fully later when they 
are needed to construct the ordered phase away from the critical point. 

The mass m of the field around the symmetric point ip = determines the conformal 
dimension of the dual operator O in the CFT 2+ i 

A± = 3/2 ± y/9/4 + m 2 (2.3) 

This also controls the asymptotic behavior of the scalar field 

^( r ) = ^A7 + ^ + --- as r ^°° ( 2 - 4 ) 
with the metric asymptotically approaching 

dv 2 

ds 2 = r 2 (-dt 2 + dxidx 1 ) + — (2.5) 



In order for us to be able to add a double trace operator as in (1.1 ) in a controlled fashion 
(without destroying the AdS^ asymptotics) we require that the mass be in the range: 

- 9/4 < m 2 < -5/4 (2.6) 

In this range, both sets of modes are normalizable and one has a choice of boundary condi- 
tions for quantizing the bulk theory. Standard quantization corresponds to setting a = 0, 
and O has dimension A+. Alternative quantization corresponds to = 0, and O has dimen- 
sion A_ [17J. We will refer to these two theories as AdS^ 6 "^ and AdS^ 1 . We want to start 
with alternative quantization, so the dimension of O is 

1/2 < A_ < 3/2 (2.7) 
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and add a double trace deformation. In this range, adding the double trace operator (1.1) 
amounts to studying the gravitational theory in asymptotically AdS± space with new bound- 
ary conditions for the scalar [151 02] 



P = na (2.8) 
Up to an overall normalization, a = (O). Note that k has dimension A + — A_ = 3 — 2A_ and 



is thus a relevant coupling in the range of interest (2.7). Hence adding the perturbation will 



induce an RG flow from the original CFT to a new theory in the IR which can be understood 



by sending k — > oo. Dividing (2.8) by k as we take this limit we see that we arrive at a = 0, 
the theory with standard boundary conditions. If we had started with K = 0we would stay 
at the unstable fixed point in alternative quantization. Many details of this flow have been 
studied (see, e.g., [3T]). 



2.2 Symmetry breaking from double trace deformations 

We will first study a simple model of spontaneous symmetry breaking. We will work with 
zero chemical potential, and hence can set the bulk Maxwell field to zero. This is natural in 



a theory which preserves charge conjugation. Adding a double trace term (1.1 ) with negative 
coefficient is expected to destabilize the vacuum. This is easily seen as follows. Consider the 
two-point function for the operator O in the vacuum (at zero chemical potential, temperature, 
and double-trace coefficient.) The retarded Green's function is (up to overall normalization) 

Gf(p)=p 2A ±- 3 p 2 = -( u + i e f+f (2.9) 

where the sign +(— ) indicates the Green's function for the standard (alternative) quantized 
theory. Also u and p are the energy and momentum respectively. If we start with the 
alternative quantized theory and add a double-trace term of the form (1.1), one findfl 



G R (p)= i l — +0(1/N 2 ). (2.10) 

T" K 



pole 



For k > 0, this just introduces a new massive pole with a width at Pp 0le = (— /t) 1 ^ 3_2A "- ) . 
However, for k < 0, we have a tachyonic instability, with a pole at real positive p 2 poie . All 
of this is directly analogous to a massless free scalar field getting a massive deformation of 
either sign. In [32] an exponentially growing tachyonic mode was explicitly found in the bulk 



2 This can be calculated either in the field theory at large N by summing up a geometric series of diagrams, 
or on the gravity side with proper treatment of boundary conditions. 
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precisely when k had the wrong sigrj^j It was widely believed that theories with k < would 
not have a stable vacuum, but in [18J it was proven that for many scalar gravity theories, 
there was a stable ground state with nonzero a = (O). 

The stability of the dual gravitational system depends on the global existence of a super- 
potential Pjjp), and the zero temperature broken symmetry ground state is given entirely 
by P c (ip)- The details can be found in [15] . The key result is the behavior of the off-shell 
potential Via). It turns out thatj^] 

V = 2(A+- A_)(W + W ), (2.11) 

W(a) is given by our boundary conditions = W'(a), and for our double trace deformation 
is simply iy(a;) = ko? /2. iyo(aO is found from a scaling limit of smooth horizonless static 
solitons in global AdS, again see [IS] for details. At T = \i = 0, scale invariance implies 



W (a) = S -^\a\^ A ' . (2.12) 

The coefficient s c depends on the full bulk potential V(ip), and is generally positive. If s c is 
negative the theory is somewhat sick since Wq is unbounded, and the alternative quantized 
theory is unstable, as it has states with arbitrarily negative energy. We do not consider this 
case any further. The full off shell potential is thus 

V(a) = (A + - A_) La 2 + 2S ^\af A A . (2.13) 



For our range of interest (2.7) the second term dominates at large a, and we have a classic 
example of spontaneous symmetry breaking with a saddle-shaped potential for negative k 
(see Fig. [T]). The ground state which minimizes (2.13) has 

/ K \A-/(3-2A_) 1 / Q_9A \ 2 

a=(o)= H) ■ v -»=-K^w) ( - k)3/ ■ < 2 - i4) 

As mentioned above, the gravitational description of this ground state is uniquely determined 
by the superpotential P c {ip)- 

Putting the theory at finite temperature can lift this instability. The detailed calculation 
is in appendix [Bj but the result is that as we heat the system up to 

T-k( withfc - 3 ( r((A + -A_)/3)r(A_/3) 2 Y /(3 - 2A - ) 

Tc -^-«) ' Wlth kl ~ 47 rr((A_ - A + )/3)r(A + /3)2 ) (2 ' 15) 



3 Strictly speaking, [32] studied the theory on a sphere, where coupling to background curvature induces a 
positive double trace term for scalars. In that case, k needed to be sufficiently negative to find the instability. 
Since we are studying the theory on Minkowski space, the critical point is simply when k changes sign. 

4 The overall factor of 2 (which was not present in [TS]) arises since we do not have a 1/2 in front of our 
action (2.1 1. 
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Figure 1: The off-shell potential as we tune k. The black dashed curve is the fine-tuned 
theory with k — 0, the blue curve is k < 0, and the red curve is k > 0. This is a strongly 
coupled version of the standard Landau- Ginzburg symmetry breaking mechanism. 

the system returns to the symmetry preserving state. Note that everything scales as a power 
of k, since this is the only scale in the problem. Another way of studying this system at 
T > is to construct the finite temperature generalization of V. This can be obtained from 
the family of hairy black holes in the bulk, as described in |33j . 

2.3 A novel holographic superconductor 

We saw above that adding a renormalizable double trace coupling can break a U{1) sym- 
metry at low temperature. This provides a new mechanism for constructing holographic 
superconductors. Unlike the previous approach, which required a nonzero charge density to 
generate a low temperature condensate, we can now work at zero net charge. In this case, 
the Maxwell field remains strictly zero, even when the charged scalar hair is present in the 
bulk. 

We computed the critical temperature above. To study the system away from T c , we 
need to specify the full nonlinear bulk potential. Working in units of k (which is analogous 
to working in units of /i or p in cases of finite density), we find that the order parameter 
behaves just as it does in the case where we find an instability by lowering Tj\i in the 
standard holographic superconductor setup. We can also calculate the difference between 
the free energy of the hairy black hole and the normal black hole, which is simply AdS- 
Schwarzschild, and find generically that below the critical temperature the hairy black hole 
is always preferred. An example which comes from a consistent string theory truncation is 
shown in Fig. [2j 

As usual, to compute the conductivity, one starts by perturbing the Maxwell field and 
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T/(-«) T /(-f) 

Figure 2: The order parameter a = (C?) and free energy density / across the second order 
phase transition down to zero temperature. In the figure on the right, the dashed red line 
is the free energy of the normal phase (Schwarzschild AdS solution) and the black line is 
the free energy of the condensed phase. We used a case with A_ = 1, and bulk potential 
V(if>) = sinh 2 (^/ v / 2)(cosh(v^V)-5) = -6-2^ 2 + C(^ 4 ) [MIES]. In [IE] it was found that 



0.56 for this potential. From (2.15), we have T c J{—k) ~ 0.62. 



metric. As shown in Appendix [A] the conductivity can be simply related to the reflection 
coefficient in a one dimensional Schrodinger problem: 

- b"(z) + V Sch (z)b(z) = cA (2.16) 

where 5A X = a x e~ tut and b = \J G(ip)a x . The Schrodinger potential is bounded (and given 



explicitly in (A.24)). z is a new radial coordinate that vanishes at infinity and goes to minus 
infinity at the horizon. To obtain the required ingoing wave boundary condition at the 
horizon, we assume b = e~ wz + TZe luJZ near z = so that b = Te~ luJZ near the horizon. The 
conductivity is simply given by 

- = \=4- (2-17) 

1 +11 ^ ' 

Since the potential is bounded, this will produce the usual behavior of the optical conductiv- 
ity. At low temperature there will typically be a gap at frequencies below the height of the 
potential and at higher frequencies the conductivity will approach its normal state value. 
There will be a delta function at u = in the condensed phase, which can be seen from a 
pole in the imaginary part of the conductivity. 

The key difference from the holographic superconductors at nonzero charge density, is 
that there is no delta function in Re [a] at u = in the normal phase. This awkward feature 
of the previous construction arose since a state with net charge can be boosted, yielding a 



5 Since the Schrodinger potential is no longer positive definite, one can sometimes get peaks at low 
frequency in the conductivity [36] . 
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nonzero current with zero applied electric field. This implies infinite DC conductivity. Since 
we can now start in a state with zero charge, we no longer have this problem. Mathematically, 
the delta function arose since the Schrodinger potential was nonzero in the normal phase 
due to a contribution from the background electric field. Here, the Schrodinger potential 
vanishes in the normal phase. The background is just the Schwarzschild AdS black hole, and 
a = 1 with no delta-function contributions. 



2.4 Non-zero density and stability conditions 

In addition to providing another way to construct holographic superconductors, the addition 
of a double trace perturbation provides a new knob for adjusting the critical temperature 
of traditional holographic superconductors. As discussed in the introduction, adding a term 
with k > makes it harder to condense the operator O and lowers the critical temperature. 
We will see below that in some cases, T c — > as k approaches a finite value, k c . This is a 
new quantum critical point which will be studied in detail in the next section. 

With a nonzero charge density, the normal phase is described by the Reissner-Nordstrom- 
AdS (RN-AdS) black hole. The critical temperature is determined by looking for a static 
normalizable mode of the scalar field in this background [13]. This marks the onset of the 
instability to form scalar hair. This problem only requires the leading terms in the functions 



V(ip), G(V'), J r (V') given in (2.2). The addition of the double trace term changes the critical 
temperature since it changes the boundary condition on the normalizable mode. 

For reference, the RN AdS black hole is described by the metric and gauge potential, 

^ 2 = -M 2 + rW + ^ * = /i-^, f = r 2 -^ + f~ 2 (2.18) 

j r 2r 4r 2 

2p 3 ,PP rp P ( 3p P 2 \ , n in , 

m o = ^r + TT> T = — — - — (2.19) 



The horizon is located at tq = p/fi, where p and p are the charge density and chemical 
potential respectively. 

On this background the linear fluctuations are given by, 

( ifi'V ( 2 2,^ 9P 2 r 2 {u + q<p) 2 \ 

[r ftp) = I m r +p - — 1 tfj (2.20) 

where we have included momentum dependence p and frequency dependence u and ^ = 
^ e -iujt+ix-p _ ^y e w jjj neec i these in Section 3 but for now they can be set to zero. Before going 
onto the case of double trace boundary conditions, we first recall some known results on 
stability conditions. One natural question to ask is what is the condition for the absence of 
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an instability at any T. In other words, when is the normal state stable at zero temperature? 
A necessary condition was identified in [T3] : the extremal RN black hole in the near horizon 



limit becomes AdS 2 x R , so if the effective mass of ip derived from (2.20) is below the AdS 2 
BF bound, the RN BH will be unstable below some critical temperature. So one condition 
for stability is demanding m 2 eS > — 1/4 where 

2 2 

9 7n q 

m eS = -^-J-9 (2.21) 

This is equivalent to demanding the conformal dimension of the operator dual to ip in the 
AdS 2 region, 5±, are real, where 



S± = \±^l + ml s . (2.22) 

This condition is too weak however, and we would like to refine it. As we will see below, 
the stability condition m^ ff > —1/4 is sufficient for standard boundary conditions for the 
scalar a = (k = oo). However since alternative boundary conditions /3 = (« = 0) 
are weaker, the scalar field can still be unstable to forming hair despite the BF bound 
in AdS 2 being satisfied. This is because for any effective mass, there are always unstable 
modes in AdS2- It is just that they are usually thrown out by the boundary conditions. With 
alternative boundary conditions in the asymptotic AdS^ region, some of these unstable modes 
are allowed. 

As further indication of the fact that alternative boundary conditions are more unstable, 
it was noticed in [T3] that T c diverges as one approaches the unitarity bound A_ = 1/2. As 
shown in Fig. 3, this divergence actually takes the form 

Tc ~ (A_ -1/2)V2 (2-23) 

Interestingly, for neutral scalar fields there is no divergence and T c approaches a finite limit 
as A_ -> 1/2. 

We now want to include the effect of nonzero k. Working at fixed [i the relevant scale 
invariant quantity that we will vary is k//i a+_a ~ (as well as T//x). It turns out to be easy to 
study T c as a function of k simply by changing the definition of "normalizable" |^] Increasing 
k always decreases T c . If the mass and charge of the scalar field is such that extreme RN AdS 
is unstable with standard boundary conditions, then T c remains nonzero for all k. However, 
if extreme RN AdS is stable with standard boundary conditions, then T c must vanish at a 



6 Amusingly this is a much simpler problem than the usual shooting problem. Rather than adjusting T 
to find a static normalizable mode we can simply fix T = T c , shoot to the boundary and read off n(T c ). 
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Figure 3: The critical temperature as a function of A_ for q = .1, .25, .5 



finite value k = k c . Both cases are illustrated in Fig. |4j Since T c can be arbitrarily large 
at k = and vanish at k c , we see that k is a very sensitive knob to adjust the critical 
temperature of the superconductor. The point Hi — Kjq is the quantum critical point that we 
will study in the next section. Notice that as k becomes large and negative in Fig. [4], fi 
becomes less important, and both curves approach the scaling (2.15). 



3 Quantum critical point 

We now turn to a more precise discussion of the point k = k c where the critical temperature 
T c /n — > 0. Below we will outline the various requirements on bulk parameters to achieve this 
critical point. Note in particular we need not fine tune these bulk parameters. The coupling 
that we do tune k/ /z a+_a - is a well defined boundary theory coupling. 



3.1 The flow of double trace couplings and the 2 point function 

It is useful to think of the extremal RN black hole as representing a flow from a CFT2+1 in 
the UV to a CFT 0+ i in the IR. This flow is induced in the UV by turning on a source (/^) for 
the charge density operator J 1 . The IR CFT can be seen by taking a scaling limit of (2.18) 
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Figure 4: The critical temperature, in units of chemical potential, as a function of the UV 
double trace coupling k for fixed A_ = 1 and q = 1/2. The top curve has g = 0.2 and has 
nonzero critical temperature for all k. The lower curve has g = —0.2 and ends at a quantum 
critical point. 



towards r — > tq at extremalityj^] Rescale coordinates as: 

f = e(r — r*) i = et x = r+x (3-1) 

where r± is the location of the horizon at extremality = r |T=o — P*l V with = p|t=o — 
/i 2 /Vl~2. Formally we can scale towards r — )• r* by expanding in e then setting e = 1. This 
yields 

ds 2 = { -6f 2 di 2 + — j + (dx 2 + df) , A = Vl2fdi (3.2) 

This is the classic AdS2 x i? 2 geometry found in the IR of many extremal black hole solutions. 
Notice in particular that the scale /i has dropped out. Since this geometry is supposed to 
be dual to a scale invariant theory, this had to be the case. The only knowledge that this 
theory has of the scale /i is encoded in higher order irrelevant terms which we have dropped. 
For example keeping the next order terms in the expansion in e one finds: 

ds 2 = ds 2 + 8 h ds\ + . . . , A = A + 8 h A x + ... (3.3) 

where ds\ and A\ have energy scaling dimension 1 under the AdS2 scaling. The chemical 
potential now appears through 5h = 1/p. Note that since 5h has dimensions of —1 this 



'This limit can be taken more carefully, keeping a finite but small T. See Appendix C. The scaling limit 
was discussed in [TT] and the discussion here follows that paper closely. 
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represents an irrelevant coupling, which when turned on induces a flow in the UV to AdS±. 
It is useful to think of 5h as opening up the R 2 directions of the metric. 

We now return to linearized fluctuations of ip in the extremal RN background. Our goal 
will be to compute the two point function of the order parameter at small frequencies and 
momenta uj,p -C fi. We proceed heuristically, leaving details to Appendix C. Following [11] 
we do a matched asymptotic expansion where we split the geometry into two regions. In 
both regions we do a perturbative expansion in e where we redefine 

lo — > eu , p — » ep (3.4) 

so as to access small frequencies and momenta. In the inner region we rescale coordinates 



as in (3.1). In the outer region we leave the coordinates unsealed. A systematic expansion 
in both regions is defined in this way, matching occurs in an intermediate region connecting 
the two. 

At zeroth order the outer region simply follows from setting u = 0, T = 0,p = in the 
full RN background. A general solution to the resulting equation can be characterized by 
the behavior at the AdS± boundary and at the extremal horizon, 

ip(r -> oo) -> a r~ A - + /3 r" A+ , ip(r ->■ r*) ->■ a {r - Tit )~ s - + /3 (r - n)' s + (3.5) 

2) -'ft) -(SO 8) (3 - 6) 

where a ± and b ± are constants (in units of /x) and can only be computed numerically^] If 



we impose linear boundary conditions (2.8) on the allowed fluctuations, then this maps into 
the following condition near r*, 

^r- = — f Q - ^—77 C ^ r where k c = b + /a + and detL = /J 2 -^ 1 ^ — (3.7) 

«o det L — a~a+(K — k c ) A + — A_ 

We would like to make the identification of /3q/&q above with the value of a double trace 
coupling in the IR AdS2 CFT, kjr. We will only be interested in k close to k c such that, 

. IR =^(.-. C ) (3.8) 

It is then natural to identify the critical point that we observed in the previous 

section. One main reason for this identification is the fact that for k < k c the double trace 
coupling in the IR is negative, and thus analogous to the discussion in Section 2, there will 
be a new state with lower free energy and scalar hair. 

8 I n [H] these same constants were called a± , b±* . 
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We are now in a position to complete the computation of the two point function of the 
order parameter. The retarded Green's function follows from imposing incoming bound- 
ary conditions at the extremal horizon in the inner region. Then to zeroth order in the e 
expansion one finds that /3q/&q = Hr(uj) where is the retarded AdS? Green's function 



for fluctuations on the background (3.2). The Green's function in the full CFT can be 
computed using any of the usual prescription [371 EE] generalized to include nonstandard 
boundary conditions. The result is up to overall normalization, 

Xr{<*>$= ° v ( \ + '^ ( ^ Z = (a+f/detL (3.9) 

where we have included higher order terms that can be important in X. These higher order 
terms always come from perturbative corrections in the outer region and are thus real. In 
contrast, £# is in general complex. The elipses above represent even higher order terms that 
we have dropped. We compute X in Appendix C. The result can be written as, 

X(uj,p) = c p p* - c u u 2 - c t k c T + c q q f-u + ^=Tq^ (3.10) 

where c, are constants in units of \x. They have the following positivity constraints depending 
on the value of <5_: c p > 0, ct > always, c w > for 5_ < —1/2 and c q > for <5_ < 0. 



The AdS2 Green's function plays the role of the self energy in (3.9) and is given by 



S R (u;,T = 0) = he^i-iu) 1 ' 25 - (3.11) 

where h is a real positive number, and e %< ^ is a phase, the precise form of which does not 
matter. We can generalize the above discussion to finite but small T. At finite temperature 
£# takes the form of a nontrivial scaling function, 

S ( uT)= (2*T\ S+ ~ S - £g± ~ *-)r(*+ - iq/V^)nS + + iq/y/3 - iu/{2ttT)) 
R{UJ ' ' V 3 J r(5^-5 + )r(6^-iq/V3)r(d- + iq/V3-iu)/(27rT)) l ' ' 

Importantly the AdS% Green's function always satisfies the constraint, 

wIm£ R (w,T)>0 (3.13) 

which is necessarily true for any bosonic spectral density. 

Generally speaking since the quantities computed in the IR AdS 2 geometry depend simply 
on two numbers q and <5L we will call these quantities "universal" . Since they are associated 
with a CFT this language seems appropriate. Other quantities that come from the outer 
region such as a* 1 , b ± and the q are "nonuniversal" , they can be computed only numerically. 
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3.2 Properties of the critical point 



Given the two point function (3.9) we can now understand the physics close to the critical 



point. First we study the phase boundary in the (k,T) plane where the order parameter 
condenses, or equivalently, where the correlation length diverges. 

We again look for a static normalizable mode at T = T c which manifests itself as a zero 



frequency pole in (3.9) . Since the instability kicks in first for the homogenous mode we 
can take p = 0. There are two cases depending on the IR conformal dimension For 
< 5- < 1/2 we can ignore X altogether, however for 5_ < the analytic correction oc T is 
larger than E# oc T 1-2 " 5 -. Thus we find for small Kir, 



(0 < 5_ < 1/2) T c 
(5-<0) T c 



k 2 



-«ir) 



1/(1-25.) 



(3.14) 
(3.15) 



where, 




r(a- - s + ) 
r(«s + -<L) 



T(5. - iq/y/Z) 



T(5 + - iq/VS) 



1/(1-25-) 



—ctk c + c q q 2 27r/ \/3 



(3.16) 



Note that while hi is a universal number (one to be compared with (2.15)) k 3 is nonuniversal 



depending on quantities ct and c q defined in the outer region. Indeed it is not clear the sign 
of &3 is fixed. Although c q > and ct > for the range of dimensions of interest, k c does 
not have a fixed sign. Generically it seems that for a critical point occurring with <5_ < 
and q = then k c < so that is fixed to be positive in such a case. However we do not 
know a proof of the positivity of ^3. 



To check the scaling relations (3.14) and (3.15), we have numerically computed the critical 



temperature. In Fig. [5] we plot T c as a function of k for various <5_. The results are perfectly 
consistent with our scaling relations. 

Now we move away from the phase boundary and study the disordered phase at zero 
temperature. Here we can examine the structure of the retarded Green's function in the 
complex u plane. We will be particularly interested in the dispersion of the mode that 
becomes tachyonic on the ordered side. Again depending on the value of 5_ the dispersion 
will differ, also now there will be a difference if the order parameter is charged or not. 
Examining (3.9) at T = for the charged case the dispersion of the pole in the complex 
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(K-K c )/fl 

Figure 5: The critical temperature close to k c for different values of <5_. These are for 
theories with q = 0, A_ = 1, and from left to right: 5_ = 0.45,0.30,0.26,0.15,0,-0.15. 
Note that when <L > the critical temperature vanishes with a power law, but for <5__ < 
it vanishes linearly. 



plane is, 



(0 < <5_ < 1/2) to* = ie 
(5- <0;q^0) L0± = UJ R (p) 



i(j>/(l-2S-) I C pP + K IR 



E R (u R (p),T = 0) 



C p p 2 + K m 



(3.17) 



(3.18) 



where in the last case the width of the quasiparticle scales like IhiEr^wr) oc IwrJ 1 " 2 ' 5 -, which 
is smaller than the energy \u R \. In this case we have a genuine quasiparticle with a mass 
oc kir. In the first case (3.17) the width scales like the energy and thus the pole does not 
represent a genuine quasiparticle. 



For the neutral case we again get the same behavior as in (3.17) however now for the 



range of conformal dimensions —1/2 < 5_ < 1/2. For the remaining neutral case, 



(<L<-l/2; 9 = 0) <4 = u R (j>f 



E R {co R (p),T = 0) 



C p p 2 + Kir 



(3.19) 



Cu 

where there are now two quasiparticles at positive and negative energies ±uj r . The width of 
these quasiparticles scales as \(jj R \~ 25 ~. 

As expected, Km. is playing the role of a mass (or an energy gap.) For all cases one can 
show via application of the constraint (3.13) that for k ir > the quasi particle pole always 
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lies in the lower half plane and there is no instability. There is a "mass gap"[^J E g , for all 
cases which is roughly given by the closest approach of the pole to to — 0. For the less 

1/2 

universal cases (3.18) and (3.19) one finds the expected results, E g oc Kir and k ir for q ^ 
and q = respectively. However for the more interesting "critical" case the gap scales as, 



E 9 ~ 



(3.20) 



The correlation length for all cases scales as £ oc '. 



When kir < 0, application of the constraint (3.13) shows that the pole always lies in the 



upper half plane for momenta p < a/— ^ir/c p , representing an instability. 

Exactly at the critical point kjr = we find a free gapless mode which disperses in the 
complex plane as 

2 

lu ~ \p\ z with z = - — (3.21) 

1 — 2o_ 

for the "critical" case, and z = 2 for g ^ and 5- < and z = 1 for q = and 5_ < 
— 1/2. We thus conclude that the physics of the critical point has a nontrivial dynamical 
critical exponent determined by the dimension of an operator 5_ in the IR AdS 2 CFT. This 
identification is consistent with the relationship E g oc £~ 2 in all cases. Interestingly z has a 
lower bound in this model, with z > 2 for the charged case and z > 1 for the neutral case. 

To summarize, in the most interesting "critical case" where < 5_ < 1/2 for q ^ and 
— 1/2 < 5_ < l/2forg = the two point function close to the critical point has the universal 
scaling form, 

XR = - , ; ^7777-77^ (3-22) 



«m + Cpp 2 + T 2 / Z g{u/T) 



where z 



1-25- 



and g(u/T) is a universal scaling function that follows from (3.12). Note 



that since the correlation function is analytic in p (the self energy is momentum independent) 
the critical point is 'locally' quantum critical [25]. This type of criticality goes beyond 
the usual Landau Ginzburg paradigm often applied to quantum critical points, due to the 
existence of the locally critical modes associated with AdS 2 . These results are compatible 
with experimental measurements of the spin susceptibility in a heavy fermion compound 
CeCu 6 _ x Au x at criticality |39j. In order to compare the spin susceptibility to the two point 
function of the triplet staggered order parameter (which is effectively Xr)i one must shift 



the momentum in (3.22) by the ordering vector associated to the anti-ferromagnetic order 
p p — K . Very similar results were also found theoretically for the spin susceptibility in 



9 This is not a genuine gap, since there will always be gapless incoherent junk coming from the AdS2 
Green's function. However it does represent a gap to the coherent part of the 2 point function, which is 
represented by the dispersing pole. 
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[25| [26] . The most important feature for comparing to experiments was u/T scaling of the 
susceptibility at the ordering vector p = K and a nontrivial exponent z ~ 2.7 or <5_ w .13. 
Indeed we capture both features here, although since our <5_ does not take a universal value, 
it is hard to make a prediction for this exponent without a real string embedding where 5- 
will be fixed. 



3.3 Renormalization group interpretation 

We would like to now give an RG interpretation of the above results and in so doing try to 
understand what to expect of the zero temperature ordered phase when k ir < 0. 

We have already discussed a major aspect of the RG flow: the extreme RN black hole 
represents a flow from one 2 + 1 CFT in the UV to a + 1 CFT in the IR. We would like to 
now understand how this picture changes in the presence of the scalar ip. For this purpose 
there is clearly a set of a distinct cases depending on the value of 5- the conformal dimension 
of the operator dual to ip in the IR CFT. Firstly for £_ complex the IR CFT will never be 
realized and there will be no critical point. We do not consider this case further here. There 
are two remaining cases < 5_ < 1/2 and <5_ < which we turn to now. 



< 5_ < 1/2 



This range of dimensions is in the "critical" range identified above where the 2 point function 



takes the more universal form (3.22) for both neutral and charged cases. Here we argue 



that this result is universally controlled by the AdS 2 theory supplemented by double trace 
deformations. The bulk field ifi in the AdS 2 x R 2 geometry is dual to a set of operators ^fp 
where p labels the momentum in the R 2 direction, which one should think of as a charge 
under the KK reduction down to AdS 2 - *&p can be viewed as the image of 0(t,p) under RG 
flow. 

The operator dimension of typ in the + 1 dimensional CFT is given by S± + 0(p 2 ), where 
in this range of dimensions we can take either value. The two point function of *f>p is then, 



(3.23) 



where the previously defined AdS 2 Green's function (3.11) is £^ = ^R- H we take the 
dimension of typ to be <5_ such that we are working in alternative quantization then we can 
reproduce the result of (3.22) at T = simply by including the following deformation of the 
+ 1 CFT [20], 

d 2 p 



5*0+1 ~~ ^ S\ 



o+i 



(3.24) 
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Note for example when the deformation is zero, Km + Cpp 2 = 0, the two point function (3.22) 
scales as u 25 '~ l consistent with typ having dimension <5_. 

At zero momentum there is a single relevant double trace coupling ktr at the critical 
point, that does not explicitly break the symmetry. Just as in the AdS^ case discussed in 
Section 2, a positive Kx& will induce a flow to standard quantization where the operators 
now have dimension 1 — 5-= 5+. On the other hand, a negative ktr will induce an instability 
which will lead to a symmetry broken state. Since c p > the zero momentum mode will go 
unstable first, so in the ordered phase the homogenous mode *$?p=o develops a vev. 

The ordered state can then be studied using gravity with the bulk field ifi turned on. For 
<5_ > 0, m? eS defined in (2.21 ) is negative which means that for the disordered phase the scalar 
ip is sitting at a maximum of an effective potential; in the ordered phase it will roll away from 
this maximum. The resulting geometry will tell us where the theory ends up in the deep 
IR. In fact there are many possibilities which will depend on the bulk functions V(ip), G{ip) 
and J(ip). One possibility which we highlight in the next section is the theory flows to a 
different AdS2 x R 2 geometry, where the field sits at the minimum of a particular effective 
potential to be discussed later. In Section 4 we will outline some other possible examples of 
deep IR geometries. At this stage to be appropriately noncommittal we call this gometry X. 

The above discussion can be understood within the context of the UV completion of 
AdS2 (a fancy name for the RN black hole) where we identify the coefficient of the double 
trace operator k as the microscopic control parameter which allows us to probe the critical 
point. The critical value k c is thus not special from the perspective of the UV theory. The 
blue region of Fig. [6] is a pictorial description of the RG flows represented by the (extremal) 
RN black hole phase of the theory. In the large- N limit, this flow is trivial, only effecting 
the fluctuations of fields in the bulk through boundary conditions. When k < K c , the double 
trace coupling runs negative in the IR theory, and the instability ensues. The full AdS± 
theory then allows us to view the end point of this instability - the geometry in the extreme 
IR flows to a new attractive fixed point X . These flows are represented by the white region 
in Fig. [6j 

Finally we make an argument as to how the order parameter behaves in the condensed 
phase. The argument 10 is based on an RG analysis close to the critical fixed point and relies 
on the picture given in Fig. [6j We will confirm the result with a numerical calculation in 
the next section. Consider shooting from the theory X in the extreme IR up to AdS±. As 
suggested by Fig. [6] there will be two tuning parameters, of which only an appropriate scale 

10 A similar argument appeared in |40j in the context of a BKT type transition, and related arguments 
appear in |41) . 
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B, x 

Figure 6: Renormalization group flow interpretation of the quantum critical point (denoted 
by the red dot) for < 5_ < 1/2. This picture is heuristic, since the couplings are only 
well defined close to each fixed point. The blue region is the disordered phase which we 
have studied in this section. A distinction is made between the theory in alternative (alt) 
and standard (std) quantizations. The white region denotes the ordered phase which we will 
study more carefully in the next section. The theories denoted by X and Y depend on bulk 
couplings. Y is the end point of the flow induced by turning on a negative double trace 
coupling in the AdS^ theory and setting fj, — 0. These flows were considered in Section 2. X 
is discussed in the text above. 



invariant ratio of the two will matter. As we tune this ratio we can shoot closer and closer 
to the AdS2 critical theory. The dashed line in Fig. [6] is an example of such a flow. In the 
AdS2 region, the metric and Maxwell field are given by (3.2) and the scalar field takes the 
form 

ip(r) ~ v f~ s - + w r~ s + (3.25) 
Here v and w are constrained by scale invariance: 



w 



(3.26) 



where s* R is a number which can be computed numerically. It is analogous to the s c in (2.12). 



We can now match (3.25) onto linear fluctuations in the extreme RN background which 
we have considered above (what we called the outer region above, see for example (3.5)). 

r = A(r - n) , t = t/A , a = vA~ s ~ , j3 = wA' 5+ (3.27) 
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where we have allowed for an arbitrary rescaling A which will drop out in the end. Addi- 
tionally we have for k close to k c , 



a = a + a , /3 = K m a (3.2£ 
Putting these results together we find, 



(O) = a » a + r -|£J " " , (3.29) 

This is precisely an IR analog of the scaling we derived in section 2 for a negative double 
trace perturbation (2.14). Note we have made an assumption that s* R is positive. This is a 



nontrivial assumption and will depend on the bulk couplings, as for example was the case 
for the analogous AdS^ parameter s c [H]. However whereas in the AdS± case, a negative s c 
meant a somewhat sick theory, we suspect a negative s* R will simply result in a first order 
transition. Since we are working with the assumption of a continuous transition we cannot 
say much about the case with s* R negative. 

More generally we can consider a black hole solution with a non-zero temperature which 
has t/> nonzero and comes close to AdS% x R 2 . Again this is a nontrivial shooting prob- 
lem. Now there is a one parameter family of solutions labeled by the temperature T and 



asymptotic to (3.25). Scale invariance imposes the following relationship between v and w, 



w = -S \-^-\ (3.30) 

where S is a scaling function with 5(0) = s* R . It can only be computed numerically. Going 
through the same matching procedure and rescaling T = AT one finds the following scaling 
relation, 

1— 24_ / rp 



~ «m = ((0)/a+) S [ m/a+)1/s . ) ( 3 - 31 ) 
This is of course consistent with the dimension of the IR operator which gets a vev 

(%) = & 5 2 (p) (3.32) 

being 5_ in alternative quantization and kir having dimensions 1 — 25-. In the next section 
we will construct these flows on the condensed side. We will confirm amongst other things 



the result (3.29) 



Actually this story is incomplete, for < 5- < 1/4 the AdS2 CFT has higher order 
multi-trace operators that are relevant (first d 4 then d 6 etc.). In this situation the scaling 
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arguments above fail, since the critical point is now multi-critical. Away from the multi- 



critical point, on the continuous side we now expect a mean fielc relationship {O) 



[— kir) 1 / 2 . The argument for this goes as follows. Nonlinear corrections 0(ip 3 ) to the linear 



equation for ip in the outer region produce additional terms in the matching (3.27). Most 
importantly 



$ = A) + Pnl + ■ ■ ■ ~ kir«o + mir«o = w & + > a = a + a NL + ■ ■ ■ ~ «o = vA - (3.33) 

where ujr can be computed along the lines of Appendix C. It is obvious that we should 
interpret uir as the RG flow of the quadruple trace operator from AdSi to AdS^ analogous 
to k — y kir. Thus at zero temperature we find, 

«mao + «m«o = s l *(a ) {1 - s - )/s - (3.34) 



For 5- > 1/4, the Uir term is not important for small do and we reproduce the result (3.29). 
However for <5_ < 1/4 the non-analytic term in a is less important and we get the mean 
field answer assuming that ujr > 0, 

(O) = a + (-k ir /u 1r ) 1/2 (3.35) 

for uir < we get a first order transition that we cannot say much about. These nonlinear 
corrections to the linear ip equation come from potential terms, as well as from back-reaction 
on gravity. They are rather complicated, however they are most likely computable in the 
probe approximation introduced in (19] for the neutral case. We leave their explicit compu- 
tation to future work. 

<L < 



In this case, defined in (2.21 ) is positive, so the IR AdS2 with ift = is stable. However, 
one can still have phase transitions which turn on ip at larger radius. In this case, we will 
see the critical exponents are not governed by the + 1 CFT but take mean field values. 



We can reproduce the more general 2 point function (3.9) for X given in (3.10) using the 
following semi-holographic action [20J, 

S = S 0+1 + J d 3 x \\d t ®\ 2 + q* + h.c) - c 2 |<9$| 2 - k$|$| 2 - ^i$|$| 4 ^ (3.36) 

+ ^/ rft (0( $ ^ + hx ) (3 ' 37) 



ii 



We thank Kristan Jensen for drawing our attention to this possibility. 
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where $ is a boundary field which we have coupled to the AdS2 CFT operator iffp and $^ is 
the spatial fourier transform of $. To reproduce (3.9) the dimension of \& must be 5 + + 0(p 2 ). 



The two point function for $ will then agree with \R with the following identifications, 



\n 



2 



1/cu , k$ = K m /c u , q® = qcqjc^ c 2 = c p /c u , $ ~ O^/cjZ (3.38) 



Note that O will also have an overlap with ^ but this will lead to a subdominant correction 



to the two point function. We have also included a nonlinear interaction term in (3.36) 
which should be generated in the flow from Ad S4 to AdS%- We will assume that > 0, but 
this need not be the case. 

In order to construct the ordered phase we first assume that we can treat rj perturbatively. 
We will show this is a consistent assumption. So for now we set rj = and work in the 
mean field approximation for $. We do this because we are working in the classical gravity 
approximation. (The whole action above should be multiplied by 1/Gn, and for small Gn 
mean field applies.) Then for k$ negative, $ develops a vev: 



($) = yJ-Kt/u* (3.39) 
Turning on 77, this will now act as a source for the homogenous mode of the operator fyp, 

s = s "» + "f¥J dt wr^ (*i + hc ) (3 - 40 > 

Since the dimension of \P is irrelevant (8 + > 1) this source will scale away in the IR. Thus 
we do not expect the vev of $ to back-react on the AdS% CFT which is now a stable fixed 
point. This is consistent with the fact that for 5- < 0, the effective mass square m 2 s > 0, 
and the bulk field ip sits at a minimum in the disordered phase and has nowhere to go in the 
ordered phase. Thus the ordering is controlled by the boundary field $. 

It is now clear how to construct the ordered phase from gravity, we simply shoot from 
AdS2 with a non-zero source term for the irrelevant operator (One also needs to turn 
on the irrelevant coupling 5h discussed around (3.3).) This situation is depicted in Fig. [7j 



Again we can match this onto the perturbations of the extreme RN black hole with 

a ~ do = vV~ K <s>/ U >5> 1 & ~ Ad + Uir&o = (3.41) 



where we have included an important nonlinear correction as in (3.33). The condition (3 = 
(which is the requirement that ip n °t blow up in the IR) and $0 = kir^o then allow us to 
match u$, 

uq,=u m Jc 2 J (3.42) 
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The vev then follows from the value of the source term do, 

(O) = a = (a + /0(-W«*) 1/2 = a + (-«m/«m) 1/2 



(3.43) 



which is same as (3.35). 



AdSf\ 



K 

AdS 4 (alt) 



AdS 2 + , 




AdS 2 + <<S»±0 



Y M- 

Figure 7: RG interpretation for 5_ < 0. In this case > 0, so from the gravity perspective 
AdS% is stable and the three fixed points on the right side are all the same. The condensation 
is controlled by an AdS% boundary field $. The three cases are distinguished by whether <3> 
has a nonzero vev (bottom), is part of the IR dynamics but does not condense (middle), or 
is massive so it drops out of the IR theory (top). 



3.4 Parametric dependence on bulk couplings 

In this section we compile some numerical results relating to the critical point. The universal 
features discussed above depend on two parameters q and 5_. However the location of the 
critical point itself k c and other nonuniversal constants appearing in the dynamic suscept- 
ability depend on three bulk parameters, m 2 ,g,q. In Fig. [8] and Fig. [9] we plot k c through 
two different slices of this three parameter space, one with q = and the other with A = 1. 



4 Constructing the ordered phase 

We have left many of the details of constructing the ordered phase to this section. We will 
consider the full back-reaction of the condensing field ip on the metric so the problem is highly 
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Figure 8: Contour plots of k c for q = 0. The lines <5_ = and k c = are shown. The solid 
(blue) region represents the excluded BF bound in AdS2 and AdS±. Note that positive k c 
tends to occur close to this region since then the theory is more unstable. For theories with 
k c negative, introducing a chemical potential fails to destabilize the theory with alternative 
boundary conditions. So these theories are more stable. 
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5- g 

Figure 9: Contour plots of k c for A_ = 1. The lines <5_ = and n c = are shown. Clearly 
increasing q tends to destabilize the theory since k c increases. Decreasing the bulk coupling 
g also tends to increase stability. 

nonlinear and we must proceed numerically. We will focus on the zero temperature case where 
the basic problem will be to find the appropriate IR solution and to integrate outwards from 
there. Perturbing around the IR solution one finds various shooting parameters that take the 
form of irrelevant couplings. These couplings generate flows in the UV to an asymptotically 
AdSi solution, representing the UV fixed point of the theory. From here we can read off 
various data such as the vev of the order parameter, the double trace couplings k and the 
free energy. 
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As one tunes the irrelevant couplings in the IR we find that one can shoot closer and 
closer to the critical point that we identified in the previous section. As we will see this 
involves a flow whose geometry is described, for a large chunk of radial proper distance, by 
the AdS2 x R 2 critical solution. 



4.1 Ansatz for the background 



We will use the following metric and field ansatz, 

dr 2 



ds z = -f(r)dt 2 + 



f(r) 



h (r)dx 



A = (p(r)dt, ip = ip(r) 



(4.1) 



This metric differs from the one used in j6] in which the radial coordinate was chosen to 
'g xx . The form (4.1) is more convenient since it allows AdS2 x R 2 as a solution. Also 



be 



numerically it is more convenient to work with this metric as we will demonstrate later. 
The equations of motion which follow from the action (2.1) are given in Appendix [A| 



(A. 5 A. 8). For our ansatz (4.1) this reduces to the system of ODEs (A. 10 A. 13). There are 



three important reparameterizations that leave the metric form invariant, thus allowing us 
to generate new solutions from old solutions. They are: 



Shift : 

Conformal rescaling : 
Spatial rescaling : 



r — 7- r + a. 

r -> Ar, (t,x) -> (t,x)/A, f -> A 2 /, 
h -> Ah, -» A0. 
x — > x/s, h — > sh. 



(4.2) 

(4.3) 
(4.4) 



We will fix these reparameterizations by demanding certain asymptotic boundary conditions, 
which we specify in the next subsection. 



4.2 Asymptotic AdS± data and the free energy 

The asymptotic UV fixed point will always be AdS±. For A + < 3 there are no irrelevant 
(single trace) operators within our truncation so AdS± is always an attractive fixed point 
in the UV. A complete expansion about AdS± can be systematically derived as discussed in 



Appendix D The leading terms include 

ip = ar~ A ~ + /3r~ A+ + . . . , = fi — pr~ l + . . . 

h = r(l + 0r~ l + ...), fh~ 2 = (1 - (m /2)r- 3 + . . .) (4.5) 

Additional terms in these expansions will be needed for some of the manipulations which 
follow, however due to their cumbersome form we leave these corrections to the appendix. 
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A general solution is parameterized by 5 constants a, fi, p, p, mo. We have used the shift 



(4.2) to fix the sub-leading r 1 terms in the metric to zero, and used the spatial rescaling in 



(4.4) to fix the normalization of h and thus the spatial components of the boundary metric. 



The remaining conformal symmetry (4.3) can be used to fix one of these five constants. 
We will work with fixed chemical potential p,, and often set p = 1. Numerical results of 
dimensionful quantities will be quoted in units of p. 

We will need to compute the thermodynamic potential, which in our ensemble will be 
the grand potential G. We will go through this in some detail by computing the Euclidean 
on-shell action. While the results for general scalar boundary conditions are known using 
other methods the inclusion of a gauge field is new. 

The grand potential is G/T = Se + S ct where Se is the Euclidean action and 5* ci are 
boundary counter terms. These counter terms are required for a good variational problem, 
as well as to regulate divergences in Se- We want to keep the metric fixed on the boundary 
(necessitating the Gibbons Hawking term), and the chemical potential fixed. For the scalar 
field ip we would like to require a, (3 to be constrained by fi = W'(a). We will work initially 
with a fixed source fi in alternative quantization where W(ai) = (3a, and generalize later. 



The counter terms in the case of fixed f3, S^, and Se are given in Appendix D with the free 
energy density evaluating to: 

gp = --m + -afiA + (A + -A_) (4.6) 

where we have defined g = G/V and V is the field theory volume. 

We can also compute how the free energy varies as we move in solution space parame- 
terized by changes in the five integration constants 5a, 6(3, Sp, Sp, 5m . Varying the on shell 
action we derive an analog of the first law of thermodynamics, 

5g p = -s5T - pdp + 2(A+ - A_)a5(3 (4.7) 



where s is the entropy density. From (4.7) it is clear that g@ is stationary at fixed T, p, and 



(3. To generalize the boundary conditions for ip we simply define 

9w = gp- 2(A + - A_)(t*/3 - W) (4.8) 

who's variation is given by, 

5 gw = -s5T - pdp - 2(A+ - A_)fo(/3 - W'(a)) (4.9) 

Notice that gw is stationary if = W(a) (and temperature and chemical potential are 
fixed). The free energy gw is what we will be concerned with in this section, the preferred 
state of the system must have lowest gw 
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To summarize, in order to have a well defined variational problem and a finite on shell ac- 
tion, for boundary conditions determined by an arbitrary function W(cx) with fixed chemical 
potential, we require the following boundary counter terms: 



= Si - 2(A + - A_) J d 3 x(a(3 ~ W) 



(4.10) 



where S^ t is given in Appendix [I5J 

At zero temperature we will not be minimizing the energy ew but rather the appropriate 
free energy at fixed chemical potential ew — PP = 9w(T = 0). To do this we present one 
final manipulation which effectively fixes the constant mo in terms of the other integration 
constants. Consider the free energy at fixed a (and T = 0): 



9a = 9/3 — 2(A + — A_)a/3 — > dg a = — pdp — 2(A + — A-)(3da 



(4.11) 



Using scale invariance we can write g a = p 3 w(a/ p A ~). Plugging this into (4.11) we find, 

w\a) = -2(A+ - A_)6 , -p/fi 2 = 3w(a) - A_aw'(a) (4.12) 

where we have defined a = a//i A " and b = (3/ fi A+ . The first equation above tells us that g a 
can be found as an integral, 



9a - 9rn = -2(A+ - A ) / da'P(a') = 2(A + - A_)W (a) 

Jo 



(4.13) 



where the integral should be understood as occurring at fixed \x (and T = 0). When we 
turn off the scalar field altogether a = 0, = there is a unique solution: the extremal RN 



black hole. This fixes the integration constant in (4.13). Using (4.8) and (4.11), the full free 
energy density is 

9w ~ 9rn = 2(A + - A_) [WoH + W{a)} (4.14) 



Notice that the right hand side is identical to the off shell potential (2.11). The above 



discussion provides a thermodynamic derivation (generalized to nonzero chemical potential) 
of this result. 

The second equation in (4.12) multiplied by p 3 becomes, 



-HP = Sg a (T = 0) + 2A_(A+ - A_)a/3 



(4.15) 



Using (4.6) and (4.11) we then find that pp = (3/2)mo- More generally with nonzero T, 
scale invariance imposes a thermodynamic relationship which can be derived analogous to 
the above manipulations. The result is, 
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mo = ^(jip + Ts) -)• g w = --(Ts + l ip) + 2(A + -AJ)\--A_aP + W(a)J (4.16) 

The first equation was derived in [12] using the Noether charge associated with the conformal 
rescaling. Practically speaking these last two equations are the most useful for reading off 
the free energy from numerical data. Since the mo coefficient is highly subleading it is more 
accurate to use p, T and s to find m using the above equation. 



4.3 IR fixed point and shooting 

To begin with we would like to construct the zero temperature ground state. One reason to 
do this is so we can see how the order parameter (O) behaves as we tune across the critical 
point. 

To understand the ground state in the ordered phase we must find the final IR fixed 
point geometry. This should not be confused with the AdSi critical point governing the 
behavior seen in section 3, but instead is the extreme infrared limit of the ordered phase we 
called X in section 3. Depending on the specifics of one's model, there are many possibilities 
for infrared fixed points, and we will not attempt to classify all possibilities here. The case 
of unbounded V{iji) was studied in [H], which found various null singularities in the zero 
temperature limit. When the bulk potential V(ip) is regulated, in the sense that it has a 
global minimum, the system is better behaved. The ground state depends on whether the 
field is charged or not. For the charged case, this problem was first considered in [12] where 
the possible IR solutions were AdS± solutions based about the new minimum or Lifshitz 
solutions. 

In this section we will mostly focus on the neutral case q — 0. Since the field ip does 
not carry any charge, there are no possible sources for the electric field lines. So there are 
one of two possible outcomes, either the horizon carries charge (the electric field is sourced 
on the horizon) or the field ip runs to a point where the coupling function G(ip) blows up 12 
Although this later case would be an interesting possibility we will leave this problem to 
future work. For now since we have a regulated potential and we are assuming that G{ip) 
is smooth then this later case will not be realized. In the former case since the horizon is 
charged, at zero temperature, there must necessarily be an AdS^ x R 2 solution in the IR, 
with constant scalar field [19]. We will refer to this solution as AdS%. 

12 Another possibility, outlined recently in [45], is that a curvature singularity develops in the IR which 
sources the gauge field, one which is physically sensible |46j . Such behavior was seen for unbounded potentials 
at q — in |44j . though not thoroughly understood. This would be an interesting state to consider since 
then the ordered phase could be insulating (gapped to charged excitations.) 
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When q = the coupling function J(ip) does not play a role. The equations of motion 
simplify greatly, reducing to 

,/ = P 

9 h 2 G(^y 



/ h J Y 2f 4//i*G(# 
h" + hi[>' 2 /2 = 0, 

h' 2 h'f hip' 2 hV{i/>) p 2 



h'f 2 2f AfWG^) ' 

where the last equation is a redundant dynamical equation. The first equation is easily 
integrated (with = at the horizon) after a solution to the other equations is obtained. A 
solution to (4.17) with constant scalar field and AdS2 x 1R 2 metric is easily found with the 
ansatz 

f = f r 2 , h = h , ^ = Vo- (4.18) 
The equations of motion reduce to 

fo = -V(iJ ), -£-(G(i/> o )Vty o )) = 0, p 2 = -2h 4 G^ )V(ij ). (4.19) 

Note that the field ip sits at a minimum of an effective potential 

Kff(V) = V(^)G(^). (4.20) 

This is our IR AdS2 solution. To flow in the UV towards AdS& we must add irrelevant 
perturbations. Since AdS% is stable, there are no relevant perturbations other than the one 
which leads to finite temperature. 

The method of studying deformations of fixed points is standard perturbation theory. 
Starting with a known fixed point $o = C0o, fo, ho), consider an expansion 

$ = $ + 5$ 1 + ( 5 2 $ 2 + ... (4.21) 

We call the perturbation relevant if the expansion is convergent at large radius and divergent 
at small radius, and irrelevant if it is divergent at large radius but convergent at small 



radius 13 Schematically, the equations of motion will be of the form 

£> 2 $i = 0, £> 2 $ 2 = $?,... (4.22) 



13 The technical definition is that if the scaling dimension of the deformation is greater than the space-time 
dimension of the CFT, it is irrelevant and classically scales away in the infrared. If it is less than the field 
theory space-time dimension it is relevant, and becomes important in the infrared. 
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One can carry out this procedure to arbitrary order to increase the precision of one's numer- 
ics. Carrying this out for AdS 2 we find 



h = ho + e{h ir ) + 0{t z 
ip = + e (ar~ A - + /3r" A + + c h x r \ + C(e 2 ), 

f = for 2 + e(a + br + c x hr 3 + c 2 ar 2 ~~ A - + c 3 /3r 2 " A+ ) + 0(e 2 ), (4.23) 



where 




A ± =^-± x r-+ " gp , (4.24) 



1 feoK//(^o) 
4 p 2 

and the constants cq, Ci, c 2 and C3 are unenlightening constants fixed by the equations of 
motion. We follow the method of identifying relevant and irrelevant deformations in [12] in 
terms of the AdS 2 scaling. Since V^A^o) > 0, A + > 1, we see perturbing the scalar field is 
irrelevant, and so we will want to turn on the mode vanishing in the infrared, namely the a 
mode. Similarly changing h is relevant, and we can turn on h x . The deformations a,b in / 
correspond to turning on a finite temperature, which we do not wish to do. With this it is 
clear that AdS 2 is a totally attractive fixed point, and we therefore expect it to be the true 
ground state. To flow to AdS^ at zero temperature we then only want to turn on a and h x , 
since they vanish in the infrared. Thanks to the conformal symmetry we in fact only have 
one shooting parameter, a/hf" , and we can use our scaling symmetries to work at fixed fi. 



Integrating (4.17) to large radius with this small perturbation at small radius we will find 



an asymptotically AdS± solution. Note that we have only presented the shooting method 



(4.23) to leading order in perturbation theory. Often (when A + 3> 1) it is necessary to work 
to higher order to have trustworthy numerical results. 

As we tune a/h x ~ we fill out a curve in the /3,a plane. We give some examples in 



Fig. 10. The basic strategy for reading off the possible states for a given W(a) is to find the 



intersection of the boundary condition curve /3 = W'(a) with the /3(a) curves. One can then 



read off the vev of the operator {O) = a and the free energy using (4.16). For example, for 



a double trace deformation W(a) = (l/2)na 2 , we look along lines = not. An example of 



the resulting (O) as a function of k is shown in Fig. 11 



4.4 Confirming the scaling relations close to the critical point 

We now examine the above results close to the critical point. This point can be seen in Fig. 
11 as the value of k = k c where (O) vanishes. As expected, this value of k c agrees well with 
the calculation of where T r = in section 2. For k > n r there are no solutions with scalar 
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Figure 10: The left figure shows (3(a) curves computed numerically as outlined above. 
The model used has V(ip) = -6 + A(A - 3)^ 2 + A^ 4 with A = 1, A = 1, and G(ip) = 
sech[(— 2gY^ 2 ip]. The curves from top to bottom represent three different choices of — g = 
.085, .125, .140. These values where chosen to give representative critical point conformal 
dimensions. The gray line shows an example double trace coupling W = na for k = —.015. 
The right figure shows the effective potential Wq + W for the same models when k = —.015. 
Note that the minima occur at the same points as the intersection points in the left figure. 
Also note that the difference in free energies in the broken phase is negative. 
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Figure 11: This figure represents a redrawing of the left panel of Fig 10 using the relation 
k = [3 /at and (O) = a. The critical point k c is where (O) — > 0. Interestingly these curves 
roughly interpolate between mean field behavior and a transition of infinite order (BKT 
transitions.) 

hair, so the ground state is simply the extremal RN black hole. For k < k c the free energy 
of the ordered phase is smaller, so this is preferred. 

First we study the metric functions of the ordered phase close to this point. The goal here 
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is to give a picture where it is clear that the bulk geometry comes close to an intermediate 
critical point. Fig. 12 shows f/r 2 and if) for flows that have k close to n c 



Three separate 

regions are clearly visible in the figure. For small r, if) is nonzero and the solution approaches 
AdS2 x R 2 . At intermediate r, if) — and the solution is approximately AdS2 x R 2 . This is 
the critical region. At large r, we approach AdS^ 



id 




ln(|A) 

-60 -40 -20 
ln(r) 



Figure 12: Metric function f/r 2 (top curves) and the scalar field if) (bottom curves) for 



different values of k close to the critical point. The model is the same as in Fig. 10 with 
a fixed g = —0.125. The two important scales are clear, the chemical potential /x which 
represents a cutoff on the critical AdS2 theory and A = (— «ir) 1-2 where 5- is the 
critical IR dimension of operator dual to if). 

We then use precision numerics to plot the order parameter and the free energy difference 



close to the critical point. These are shown in Fig. 13 for different values of the scaling 



dimension 5- given in (2.22). This scaling dimension is distinct from A_, and its importance 



was made clear in Section 3. As Fig. 13 shows, the data fits well to the following scaling 

relations to within 5% error (this error could be reduced by going closer to the critical point): 

<?- i 
(O) oc (k c - k) 1 " 25 - , g w - g RN oc (« c - «) 1 ~ 2S - 



(4.25) 



The former was predicted in (3.29), and the latter can be derived from the former using 



(4.13) 



Finally one point we would like to emphasize is that at the end of the day the details 
of bulk coupling functions V, G, J and the IR geometry X will not be important for un- 
derstanding physics near the critical point. Although, as already mentioned, there are a 



14 Similar "3 shelf" structure was seen in |47j for a probe brane model. It would be interesting to understand 
better the critical theory involved in that case. 
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Figure 13: The left figure shows the order parameter close to the critical point for the three 



examples considered previously in Fig. 10 (with the order of the curves flipped.) For large 
negative k the data all follow the same curve determined by the UV conformal dimension 
A_ . However for k — > n c the curves diverge depending on their critical conformal dimension, 
<5_ = .26, .37, .45. The right figure shows the free energy difference. 

few consistency conditions which depend on the bulk coupling functions and are required 
in order to realize this continuous second order transition. These are the conditions s c > 
defined in (2.12) and s* R > defined in (3.26). Also when it is important we require Mir > 



defined in (3.33). Although we have not attempted to systematically compute these, as we 



have shown above it is not hard to find models in which these conditions are satisfied. 
4.5 Critical solution 

The size of the intermediate AdS2 region grows without bound clS 7 f^ir 

from below. At the 

critical point, the asymptotic AdS± part of the solution completely decouples from the AdS 2 
IR region. There are two limiting solutions: One is the standard extreme RN AdS black hole 
which keeps the asymptotic AdS^ region, and the other is a solution which asymptotically 
approaches AdS2 x R 2 in the UV. 

Since we can compactify R 2 into T 2 , this seems to contradict earlier arguments that there 
are no nontrivial solutions with asymptotically AdS2 (times compact) boundary conditions 
One way to see the potential problem is to consider a congruence of null geodesies at 



constant x in the metric (4.1 ). Their convergence is c oc ±h' /h. As long as h is not constant 
everywhere, either the left or right moving geodesies will have c > somewhere. However, 
the Raychaudhuri equation and the null energy condition imply that if c > at one point, 
it must diverge in finite affine parameter causing a spacetime singularity. The resolution is 
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that our limiting solution does have a singularity, but it is inside the Poincare horizon of the 
AdS 2 IR region. In fact, its Penrose diagram is similar to the extreme RN AdS solution. 
In short, the obstruction [H] applies only to solutions which are asymptotically AdS2 in 
both the left and right asymptotic regions, and our critical solution has only one asymptotic 
region. 

We have focussed on a neutral scalar field in this section for simplicity. If q ^ 0, then 
the natural IR geometry is a new AdS 4 geometry with ip sitting at the global minimum of 
the potential. In this case, the critical solution flows from AdS2 in the UV to AdS 4 in the 
IR, just the opposite of the standard RN AdS case. 



5 Discussion 

In the growing literature on applications of gauge/gravity duality to condensed matter, one 
class of relevant operators has been largely ignored: multi-trace deformations. We have 
started to remedy the situation by studying various effects of these deformations. This in- 
cludes spontaneous symmetry breaking and new quantum critical points. Below we first 
summarize the behavior near the quantum critical points, and then discuss various applica- 
tions and generalizations of our results. 



5.1 Summary of critical exponents 



We have derived various scaling relations which hold near the quantum critical point which 
we now summarize. Let the operator dual to our bulk scalar, O, have dimension A_, and 
k be the coefficient of a double trace perturbation as in (1.1). Then for \x = 0, the critical 
point is k 



and near this point (from (2.14) and (2.15)): 

(O) oc (- K ) A -/( 3 ~ 2A -) , T c oc 



_ xl/(3-2A_) 



(5.1) 



For fi 7^ 0, the critical point is typically at a nonzero k c . For k close to k c the bulk 
solution has a large intermediate AdS2 x R 2 region. Let <5_ be the dimension of the operator 
dual to the bulk scalar in the corresponding + 1 dimensional CFT. There are a few cases 
depending on S-. For 1/4 < 5_ < 1/2, the critical point can be viewed as turning on a 
double trace operator with negative coefficient in the + 1 CFT. The exponents are given 



by (3.29) and (3.14) which are just the IR analogs of the ones above. Setting k jr oc k — k c : 

(O) oc (- KlR ) 5 -/( 1 - 25 -) , T c oc (- KlR ) 1 /( 1 - 25 -) (5.2) 
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For < 5- < 1/4, there are relevant higher multi-trace deformations. If the phase transition 
remains second order, then (3.35) gives 



(O) oc (-k 1r ) 1/2 , T c oc (- KlR ) 1 /d-^) (5 _ 3) 



For <5_ < 0, we have mean field behavior (from (3.43) and (3.15)) 

(O) oc (-kir) 1/2 , T c oc (-Km) (5.4) 

Note that the exponent of (O) is continuous at 5 = 1/4, and the exponent of T c is continuous 
at 5^ = 0. 

We have also found that at the critical point, Km = 0, there is a free gap less mode which 
satisfies u ~ \p \ z where the dynamical critical exponent z is given by 

z = 2 q ^ and 5_ < (5.5) 

z = l q = and 5_ < -1/2 (5.6) 

2 

2 = — — otherwise (5.7) 

1 — 2o_ 

5.2 Applications of our results 

As we discussed in section two, one main application of our results (with q ^ 0) is to holo- 
graphic superconductors. Since double trace deformations can break a £7(1) symmetry even 
with zero net charge density, they provide a new mechanism for constructing gravitational 
duals of superconductors. One appealing aspect of the new construction is that, despite 
translation invariance, the DC conductivity in the normal phase is finite. We also showed 
that adding a double trace deformation to the traditional construction of holographic super- 
conductors introduces a sensitive knob for adjusting the critical temperature. 

Another possible application is to neutral order parameters (neutral under the U(l) 
charge density). As mentioned in the introduction these can be used to model the onset 
of anti-ferromagnetic order [TS]. Interestingly we found examples where we could drive the 
critical temperature to zero, revealing a quantum critical point characterized by the disap- 
pearance of anti-ferromagnetic order. The nature of the critical point is rather mysterious, 



however it's locally quantum nature discussed around (3.22) has appeared in previous theo- 
retical studies of heavy fermion criticality j25j [26] . Indeed measurements |39j of the dynamic 
susceptibility of a certain heavy fermion material CeCu§_ x Au x close to criticality seem con- 
sistent with the form of the two point function that we derive. Optimistically our results 
based on AdS^ x R 2 might shed some light on the robustness of the large dimension limit used 



38 



to justify the theoretical results of [25] 126] . Or turning this the other way, we might learn 
something about the mysterious CFT dual to AdS 2 x R 2 . Such possibilities were discussed 
recently in a related context [19] . The relationship of AdS 2 to local quantum criticality was 
discussed in [50] . 

We now would like to discuss some immediate generalizations. 



5.3 Magnetic fields 



Studying the theory when a magnetic field is turned on is often difficult, due to the formation 
of superconducting droplets at nonzero q. However, we can study the theory near the critical 
point following [1, 5T] by studying linearized analysis around a dyonic black hole. The dyonic 
RN AdS black hokf^j is a simple generalization of (2.18), which now we are writing in terms 
of the horizon radius tq, the chemical potential /x = p/tq, and the magnetization M = B/ro 
(not to confused with the mass of the scalar field m) 



f 



ds 2 



2 m , 

r - 27 + 



-fdt 2 + r 2 (dx 2 + dy 2 ) + 
(/i 2 + M 2 )r 2 



m 



dr 2 

T 

(/i 2 + M 2 



A = ti(l 



1 
r 



T 



dt + Mroxdy, (5. 
/i 2 + M 2 



3r 

47T 



(5.9) 



4r 2 1 "~ v 2 ' u ' 4tt 167rr 

We can generalize our calculations of the critical temperature by studying linearized fluctu- 
ations, a simple generalization of [U EU IT9] . 

For now let us restrict to q = 0, and study the effect of a nontrivial G. The important 
equation is the generalization of (2.20), with ^ = R(r)e tky ~ tuJt , which leads to the wave 
equation 

g(v 2 



(r 2 fR' 



[roj) 



+ 



M 2 )r 2 



m r 



R = 0. 



(5.10) 



/ ' 2r 2 

This can be solved either numerically or by a matched asymptotic expansion. Since we 
are working at q = 0, turning on a magnetic field effectively shifts g for our linearized 
fluctuations. If we restrict to zero temperature and frequency, we find that the effective 
AdS 2 mass m e ff changes as a function of M, 



rn 



m 2 f M 2 - fi 2 
If + 9 \M 2 + n 2 



(5.11) 



where again S± = \ ± \l\ +m 2 ff is the infrared dimension. It is clear that turning on a 



"eff 

magnetic field changes the infrared dimension of the operator, and can possibly change it 
15 See [55] for a detailed study of the transport properties of the dyonic black hole. 
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from real to complex 16 This defines a critical magnetization where m 2 



eff 



-1/4, 



Mr 



/M 



12g - 2m 2 - 3 



(5.12) 

12g + 2m 2 + 3 V ; 

where there is an infinite order holographic BKT transition^ just as that found in [4"0l 119] 
and studied further in [53l|5ll|55] (for a field theoretic discussion see |56j.) The critical point 
at k = k c turns into a second order phase boundary, which terminates at M c . See figure 14 
for example phase diagrams. 
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Figure 14: Phase diagrams when turning on a magnetic field. On the left, when g is negative 
and q is zero, turning on a magnetic field helps destabilize the disordered phase. On the right, 
when q is nonzero and g is zero, turning on a magnetic field helps stabilize the disordered 
phase. The red dashed line is the holgraphic BKT transition at M — M c , an infinite order 
transition. The solid blue line is the non-mean field transition controlled by Km changing 
sign, where 5_ < 1/2. At the holographic BKT line the critical exponents diverge. 



Along our k c transition line, moving towards the BKT transition, <5_ — > 1/2 which implies 



the critical exponents of the transition (5.2, 5.7) diverge. This must match onto the known 
holographic BKT behavior [TjJ], where 



T c ,(0) ~exp [-Tt/^3/2 



6m 2 eff 



(5.13) 



It would be interesting to understand the behavior of the order parameter two point func- 
tion as we approach the holographic BKT transition, and how it matches on to the known 



behavior along the second order transition line. At nonzero q, (5.11) generalizes to 

,2 



m 



eff 



m 
~6~ 



+ 



qM 



3(M 2 + fi 2 ) ^3(M 2 + n 2 ) 



9 



M 2 -fi 2 
M 2 + /j, 2 



(5.14) 



16 For q this was pointed out in [E5], and we see that a nontrivial G can cause this to occur for q = 
as well. The magnetic field at which this occurs when G = I, q ^ is the B c found in [B], 

17 For a different type of quantum phase transition that occurs as a function of a magnetic field, see |S5] , 
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5.4 Quantum corrections 



We will now discuss the effect of quantum corrections governed by an expansion in the 
gravitational coupling Gn- These correspond to 1/N corrections in the field theory. We will 
be particularly interested in how robust our results for the critical exponents are to bulk 
quantum corrections. 

We start by considering the disordered phase. Firstly it would be interesting to under- 
stand one loop corrections to the thermodynamics and transport that come from fluctuations 
of ip. Methods developed in J57J [53 [59J [601 (El EE] would be useful for this purpose. In the 
large N limit, quite often bulk geometric contributions to thermodynamics and transport 
dominate since they are of order 1/Gjv, so 1 loop contributions which are of order 1 in the 
Gn expansion will be suppressed. However the mode that becomes gapless at the critical 
point will give large contributions that may be isolated from the geometric contributions 
due to strong IR behavior. This philosophy was used in [571 ES EH] and in [12] to extract 
thermodynamics (most notably quantum oscillations) and charge transport of holographic 
non-Fermi liquids. 

It is also important to understand how these modes renormalize bulk couplings when 
running in loops. Although this would require a fairly complicated bulk one loop calculation, 
we can argue what the outcome should be, based on the low energy semi-holographic model 
that we wrote down in (3.36). Note that in the disordered phase this model reproduces the 
2 point function \R f° r an <5- re& l- Consider here the neutral case q = 0, then ignoring the 
coupling 1] between the boundary mode $ and the CFT operator \l/ what we wrote done in 
(3.36 ) was simply a $ 4 theory in 2 + 1 dimensions. Thus we expect large quantum corrections 
to mean field at low energies, since m$ should run to strong coupling (at the Wilson Fisher 
fixed point). Low energies should be u ~ G^ in this case. 

We believe the real situation when r\ is non-zero will be slightly different since the modes 
in the 2 point function for $ now have a different critical exponent z = max(l, 2/(1 — 25_)) 



is 



This changes the RG power counting of the non-linear coupling $ 4 effectively softening the 
running. One can express this in terms of an effective dimensionality d e g = 2 + z, where for 
d e g > 4 one is above the upper critical dimension of $ 4 theory and non-linear effects should 
not be important (except in the case that they are dangerously irrelevant, which here they 
are not.) Thus the scaling exponents we derived should not receive corrections as long as 
z > 2. Interestingly when z < 2 we should be able to see a Wilson Fisher type fixed point 
as an expansion in 5_, analogous to the classic d — 4 expansion since d e g ~ 4 + 45_ for small 
<5_ < 0. This will then change the mean field predictions we found above, for example in 

18 This discussion was inspired by similar considerations in [25] , 
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(3.43). We leave details for future investigation. This softening of nonlinear interactions 
occurs in a similar context for holographic non-fermi liquids, leading to suppression of the 
BCS instability [60]. 

On the ordered side, of course the second order mean field transition that occurs along 
the line T = T c (k) will not survive corrections beyond mean field theory. In 2 + 1 dimensions 
at finite temperature fluctuations of a continuous order parameter destroy long range order. 
This is a consequence of the Coleman-Mermin- Wagner theorem studied in a holographic 
context in [51] . Ultimately on very long distance scales (~ c*^ Gn ) this mean field transition 
should be replaced by a (conventional) BKT transition. It would be very interesting to 
understand how these corrections scale close to the critical point. At zero temperature since 
fluctuations in the temporal direction become important it is now possible to have continuous 
order. (One should not be fooled by the fact that the critical theory seems to be controlled 
by a + 1 dimensional theory. The IR of the zero temperature ordered phase has nothing to 
do with this + 1 dimensional theory, and besides from the bulk perspective the R 2 directions 
of AdS2 x R 2 are infinite and still allow for fluctuations.) 

5.5 Lifshitz normal phase 

One possible generalization of the critical point story is to replace the normal phase, which for 
us was the RN black hole, with a more general charged black hole solution, which is however 
still asymptotically AdS^. There are several motivations for doing this, one is that this may 
solve the finite entropy problem at extremality, that plagues the RN black hole. In this vein 
one expects quantum corrections (orthogonal to the ones discussed above) to lift the finite 
degeneracy associated with this finite entropy. In so doing, these corrections must replace 
the AdS2 x R 2 critical solution with something else. This is exactly what happens when 
one studies the back reaction of fermions in the bulk [62J, where the IR AdS? x R 2 solution 
is replaced by a Lifshitz geometry [63] with a dynamical critical exponent zq ~ 0(1/Gn) 
(we use a subscript G to denote a geometric exponent, to be contrasted with our z.) For 
zq — y oo one recovers the AdS2 x R 2 geometry. So finite zq is playing the role of a regulator 
and we expect this to be a general phenomenon. Note for example now the entropy scales 
as S ~ T 2 I Z ° at low temperatures .611 . 

Another motivation is simply to have a more general understanding of what holographic 
models can do at finite density. One simple generalization is to add a dilaton field a which 
couples to the gauge field, considered first in the AdS/CMT context in [64], 19 We will 

?2 v. „2aa rp2 



consider this model here as a concrete example. In our action (2.1 ) one replaces F — > e F 



19 



See [S] for further generalizations. 
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and adds a kinetic term for a with normalization — 2(Vcr) 2 . We will not add a potential for 
a. There are many choices for how to couple o to ip, we begin for simplicity by assuming no 
extra couplings between a and ip other than through the replacement above. However we 
will argue that other choices of couplings should produce similar results. 

The normal phase has ip = 0. The appropriate solution which is charged and asymp- 
totically AdS^ in the UV was identified in [M]- It has a Lifshitz like geometry in the IR, 
replacing the AdS2 x R 2 solution of the RN black hole. Roughly speaking the dilaton runs 
a — > oo such that the coupling and the charge density redshift away in the IR. So there 
is no need for a charge on the horizon to source the electric field. The critical exponent 
is determined by the dilatonic coupling zq = 1 + 2/a 2 . The IR scaling geometry at finite 
temperature is, 

ds 2 = -fdt 2 + ^J + f2/ZGrff2 ' / = A/r 2 (l - {f T /f) 1+2 ' ZG ) (5.15) 

a = y z G~ l log(r///) , A = Mfi- 2/zG (r 1+2 ' ZG - r l T +2/za ) dt (5.16) 
z G \ J 

where J\f,A4 are numbers determined by z and fj- = (2/3)7r(l + 1/zg)T. This solution 



should be compared to (C.2). Note that the solution is not precisely scale invariant due to 
the form of both the gauge field and the dilaton, to emphasize this we have included explicit 
factors of \x. 

In this background we can now repeat our construction of the critical point. We will only 
consider the disordered phase, where the 2 point function determines most of the physics. 



The linear fluctuations of ip on (5.15) have a scale invariant form for r <C fi, with dimensions 
[cj] = zc[p\ = [T] = [r]. This is a nontrivial statement because the background is not scale 
invariant, so in the limit f <C fi various couplings disappear from the equation for ip. For 
example since the gauge potential qcj) redshifts away compared to the energy uj, the coupling 
q will not play a role in IR physics. Including more general couplings between a and ip will 
result in different IR scaling limits; the only assumption one needs in order to extend these 
results, is that this limit exists (and is nontrivial.)^ 



Imposing incoming boundary conditions on the solution to this scale invariant differential 



equation, for large f one finds a generalization of (C.3). 

V>(f) ~ (f- 5 - + Z R {u,T,p)) = (r- s - +T 5 +- 5 -g (u/T,p 2 /T 2 / ZG ) f~ 5 +) (5.17) 
where E# is the IR retarded Green's function, which now has general nonanalytic dependence 



One may also need to generalize the scaling to [uj] — zq[p\ = n[T] = n[f] for some 
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on p 2 . The critical conformal dimensions are given by, 



>± 



L.2 , f d eA 2 ^2 d% ( A(A - 3) / , 1\ /, 1 



where df s = 2 + zq. Note in particular any dependence on q has vanished. We do not 
however expect the expression for m e g to be universal when considering other couplings of 
if) to the dilaton. Matching to the outer region we find for the IR limit of the two point 
function, 

Z 

XR ~ «m + erf? - T*/*g(u>/T,iP/T*/»o) (5 ' 19) 

where z = 2/(5+ — <5_). Note that this result can also be obtained semi-holographically 
|20j . Again there should be some matching between Kjr and (k — k c ) in the outer region. 
Importantly we have included an analytic in p 2 correction which will also come from the 
outer region. This correction will be important if z < zq as we discuss now. 

The coherent part of the two point function comes from the dispersing mode with p 2 ~ 
T 2 / z ~ uj 2 / z . Along this mode the second argument in the scaling function g vanishes, 
fi/T 2 ' z G -> 0. That is, we can ignore the momentum dependence in the scaling function 
g. Thus we are forced back to a locally quantum critical theory, where the 2 point function 
is effectively analytic in p. The critical exponent z is still determined by the conformal 
dimension of the operator dual to if> in the Lifshitz IR theory. 

Of course this analyticity in momentum and the critical exponent z will not persist away 
from the coherent peak, and the 2 point function will have signatures of the z G Lifshitz 
scaling for p 2 ~ uj 2 I Zg . These signatures will be incoherent and subdominant. They are 
analogous to the gapless modes that are present in the AdS-z case, and also present in the 
holographic no n- Fermi liquids [llj. This makes sense in the context of the limit z G — > oo 
where one recovers AdS2- 

For z > z G then the analytic p 2 term is subdominant and the two point function will 
have strong nonanalyticity in both frequency and momentum. The critical point is then 
truly governed by the Lifshitz geometry. 

So to summarize, for the case z < zq we see that our results for the AdS2 critical point 
are robust under generalization to a Lifshitz IR scaling geometry. Further it is easy to see 
that the phase boundary is still controlled by the conformal dimension T c ~ (— Kir) 1 '' <5+-<5_ ^ 
and the order parameter turns on as (O) ~ (— Kir) <5_ ^ <5+_5_ ' ) which is similar to what we 
found earlier, although now we have 5 + + 5_ = d^ s /zc- A more thorough discussion is in 
order, but we leave this to future work. 
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A Equations of motion 

A.I Field equations 

We begin with the action 

S = J d'x^g (r - \g^)F 2 - (VV) 2 - JW(V6 - qAf - Vty)) , (A.I) 

where we have decomposed a complex scalar \1> = ipe 16 . We have normalized A such that 

G=l + # 2 + 0(V> 4 ), (A.2) 
and periodicity of 6 ~ 6 + 2n implies 

J = ^ 2 + 0{^). (A.3) 

Lastly we assume 

l/ = -6 + mV + 0(^), (A.4) 
where we have chosen units where the AdS radius is unity. The wave equations for the scalar 



are 



V 2 ^ - \YW) - \g'^)F 2 - \j'mV0 - qAf = 0, (A.5) 

V,[JW(V"e-?#)]=0. (A.6) 

Maxwell's equation is 

V M [G(V>)F^] + 2qJ(ij)(y u 9 - qA v ) = 0. (A.7) 

Einstein's equation is 

R,u - \g,u (r - \g^)f 2 - (WO 2 - J(V)(V0 - qAf - vy,)) 

—GWFwFS - V^V.V - J(^)(V^ - qA^)(V„e - qA v ) = 0. (A.8) 
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A. 2 Our ansatz 



When studying static black hole solutions, we make the ansatz 



ds 2 = -f{r)dt 2 + dr 2 /f{r) + h{rf{dx 2 + dy 2 ), 

A = <j)(r)dt, ip = ip(r), 6 = 0. (A.9) 

This leads to the equations of motion 

^ + g + £),_m + <^ + ^=o, (A ,o) 

* + It + -gw"J*- g(w* = 0, (A ' n) 

ft . + £ fc+ ^ k=01 (A , 2) 

fc" /w fry/ 2 i i ?M#;_ n , A1 _. 

~h f 2~ + + '~f 2f*~~ " ( j 
There is a third component of Einstein's equations which is a dynamical equation for / and 
can be derived from the above equations, 

r + ^ + yw _g(M!-^W = o. (A . 14) 

The equations of motion are invariant under the following symmetries, which leave invariant 
the scalar field, the line element, and the Maxwell one-form: First, a shift symmetry, 

r — > r + a. (A. 15) 

Second, a conformal rescaling, 

r ->■ Ar, (t,Xi) -> (t,Xi)/A, f ->■ A 2 /, h ^ Ah, ->■ A0. (A. 16) 
Lastly, a spacial rescaling 

Xi—>Xi/s, h —> sh. (A. 17) 
A. 3 Calculating conductivity 

To calculate conductivity in this background, we study linearized modes 

5A X = a x (r)e-^, 6g tx = g tx (r)e~ wt . (A.18) 
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These have the following equations of motion 



a x + 



f + G&) ) a ^ 
2ti 

9tx - ~T9tx + G(ip)(p'a x = 0. 



io 2 _ 2q 2 J(^) 



<*x + J9 ta 



2h'. 



fh 



~9ta 



We can eliminate g tx and find 



f Gty) J 



a' + 



u?_ _ G{^j)(f)' 2 _ 2q 2 J(ijj) \ 

P 7 JgW) 



0. 



(A.19) 



(A.20) 



There is now a factor of G(ip) multiplying the Maxwell term in the action, but since ip = 
a/r A - + . . . and A_ > 1/2, G goes to one rapidly enough at the boundary that the analysis 
of [6] follows and we can read off the conductivity from the asymptotics of a x , 



4°)+aW/r + 



(i) 



a 



—ia^'/uja 



(0) 



(A.21) 



As in jH] we can rewrite this as Schrodinger's equation. Using a new radial variable and 
rescaling the gauge field 



dz 



dr 

7' 



VgWY 



we find that (A.20) reduces to 



b"(z) + V Sch (z)b(z) = cu\ 



where 



V(z) = f G(^W r + 



2q 2 J(i)) 1 



2^/GW) 



(A.22) 



(A.23) 



(A.24) 



We can relate conductivity to the reflection amplitude. In the new coordinates, near the 
boundary at z — 0, where z = — 1/r + . . ., we assume h = e~ luJZ + TZe luJZ and near the horizon 
at z = — oo ingoing boundary conditions require b = Te~ luJz . We find 



a 



1-11 

l + W 



(A.25) 



We once again see that the necessary condition to find a hard gap in the real conductivity 
is that V(z) not vanish on the horizon. In jH] it was proven that for J = ip 2 ,G = 1 the 
Schrodinger potential must vanish on the horizon. It is an interesting question as to whether 
one can evade this theorem with more general J and G. 
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B Critical temperature at zero chemical potential 



We can find the critical temperature by looking for static linearized scalar hair on AdS 
satisfying our linear boundary conditions. In the background of a d + 1 dimensional AdS- 
Schwarzschild black hole, the wave equation is 

,„ //' d-l\ „ A(A-d) , (J 1 k 2 \ , N 

where / = r 2 [l — (r /r) d ] and T = rfr /47r. We have dropped the subscript and are working 
with A = A_. The static case k = u = can be solved via hypergeometric functions. 
Writing z = r /r, 

$ = Cl z A 2 F 1 [A/d,A/d,2A/d;z d }+c 2 (A <*d-A) (B.2) 
which at large r falls off as 

^ = c 1 (r /r) A + c 2 (r /r) d - A + ... (B.3) 
and therefore satisfies the linear boundary condition 

^ = K =^(47rT/d) d - 2A (B.4) 
a ci 

Now recalling that 2^1 [a, b, a + b; y] oc log(l — y) we must choose the appropriate linear 
combination to cancel the logarithmic term for regularity on the horizon (this logarithmic 
branch is coming from the u = limit of ingoing and outgoing waves on the horizon.) The 
appropriate combination is 

C2 r(2A/rf)r(i-A/^) 2 

ci " T(2 - 2A/d)T(A/d) 2 { ' 

This tells us that there is a second order phase transition as we heat the system up, and 
above the critical temperature 

lc ~4n{ T(l-A/dyT(-l + 2A/d)J [ K) ^ 

the system returns to the symmetry preserving state. 
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C More on the 2 point function 



Here we give details of the computation of the two point function in the ordered phase that 
we left out from Section 3. Most importantly we include a discussion of finite but small 
temperature, and compute the quantity X given in (3.10). 

The small frequency limit can be supplemented with a small T limit, where u/T is held 
fixed. To do this, as in (3.4), we redefine, 

oj —> eu, T — > eT, p — >■ ep (CI) 

and proceed to expand in e (setting e = 1 at the end.) The inner region now has the metric 
of the AdS 2 black hole: 

dsl = -fdi 2 + -^ + dx 2 , f = 6f(f - (2/3)ttT) , A = VV2(r - (2/3)nT)di (C.2) 

In this background we can compute the AdS2 Green's function for if), imposing incoming 
boundary conditions at the horizon. Asymptotically at the AdS 2 boundary one has, 

Mr) = M (r 5 - + E R {u, T)f- 5 +) (C.3) 



The answer is found in Appendix D Eq. (27) of [TT]. It is reproduced here in (3.12) 
using our current conventions. Higher order corrections to the inner region will never be 
important, so we do not consider them. 

In the outer region we expand if)(r) = ipo(r) + eif)\(r) + e 2 if) 2 (r) + • • •• Where: 

Dif) = 0, LhJ> 1 = X l (i/; ), Lhp2 = X 2 (Mi) ... (C.4) 
Where at zeroth order we have, 

D = d r (r 2 f d r )-m 2 r 2 + 9 ^ + q 2 r 2 (f) 2 /f , <h = f*-^, fo = r 2 - ^ + || (C.5) 
where we recall that po = /i 2 /v^l2- And for the next two orders the source terms are: 



X^o) = -~D(rif)' ) + 2q f- U + ^=Tq\ (C.6) 



X 2 (iP , = (V ~ + C(u, T)ipi (C.7) 

where in X 2 we have ignored nonleading corrections. These corrections will be of order 
qu 2 , qcoT, T 2 , which are sub- leading compared to X 1 . Note that we have kept an u 2 term in 
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X 2 since for the neutral case this is the leading correction in cu. Given this we can now treat 
ipi + eip2 ~ i^l together (L stands for leading corrections). So we now write, 



D 



2ttT 
Vo/i 



X L i> , X L = 2q 



-u 



2tt 

7! 



Tq 



fo 



+ e 



f 



u?r 2 



fo 



The horizon and boundary behavior of ipQ was given in (3.5). Since we can always add a 
homogenous solution to ipi we demand that at the horizon ipi has the following behavior 



1>i 



B 2 (r - r*)- tf+-2 + B x (r - n)~ 5+ ~ l + ■ (r - r*)~ 5 + 
A 2 [r - n)- 5 - 2 + A x {r - n)- 6 — 1 + • (r - n)' 5 - +... 



(C.9) 
(CIO) 



where A lj2 and -B 1>2 are completely determined by ?/>o through the source terms in (C.4). The 
nontrivial requirement relates to the terms we have set to zero above. Given this constraint 
we can then match to the inner region via (C.3) and (3.5), we find Pq/q.q = ^r with no 
contribution from ipL (and here we have set e = 1.) Note the diverging terms in (C.9) can 
be matched to lower order 1/f corrections to (C.3). 
Asymptotically at the boundary of AdS^ we have, 



if>L = oi L r A ~ + f3 L r A ^ 



and from this and (3.5) we can read off the 2 point function: 

a + tdi + • • • 



a 



XR 



(C.11) 



(C.12) 



-(3 + Ka —p + Ka + e(—(3 L + Ka L ) + . . . 

Close to the critical point where k ~ k c the leading correction comes only from — Pl+KcCXl 
which we intend to compute now. Consider the normalizable mode ip Kc at the critical point 
defined as: 



D^ Kc = , ij) KB = l(T- r )- 5 - + 0(r - r o y 5 + , 



a (r 



k c .t 



(C.13) 



where the last two equations are the boundary conditions at the horizon and AdS$ boundary 
respectively. Then if we multiply the first equation in (C.8) by ip Kc and integrate by parts 
twice we find, 



W 



27rT ,/ 
Vo/i 



(C.14) 



where W[ipA'i iPb] — fo^i^A^'s ~ ^b^'a) i s the Wronskian. Also close to the critical point 
Yj R is small so for the purposes of computing this correction we can set p Q = such that 
■00 = Ao^Kc- Evaluating the Wronskian we find the desired result, 



-p L + K c a L 
d 



A>=0 



2ttT 



a + (A+ - A_)k c + 



1 

a+(A + - A_ 



dr^ Kc X L *l> Kc (C.15) 



r*(reg) 
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where the integral is understood to be regulated. Possible divergences come only from the 



horizon and the regulators actually come from the Wronskian in (C.14) evaluated as r — > r 



when one includes the diverging terms given in (C.9). 
Finally the leading contribution is, 



-fa + Ka 
a 



cl + (k - K c ) 



detL. 



(C.16) 



Compiling all of these results and using (C.12) we arrive at the result quoted in (3.9) where 
the constants in X are given by, 



2vr 



det L y/3 fj, 

(g xx ), c 



(A- 



A. 



I»*l> 



where the angle brackets define the following: 

(V) 



detL(A H 



A_ 



[A t \g tt \) 



-gi\) Kc Yi\) Kc 



(C.17) 
(C.18) 



(C.19) 



r*(reg) 



and where the metric components above are all for the extremal black hole metric. Note 
that these integrals are not always divergent, depending on the value of <5_. When they are 



not divergent since all the integrands are positive it follows that the constants in (C.17) are 



positive. This is always the case for c p , it is only the case for c w when <5__ < —1/2 and it is 
only the case for c q when <5_ < 0. 



D Complete AdS^ expansion and boundary terms 



We give here the asymptotic solution to (A. 10 - A. 13) with V = —6 + m 2 tp 2 + (\/2)ip A . The 
expansion will be complete up to and including C(r~ 3 ) assuming that 3/2 < A + < 3. To 
this order, it suffices to take G = 1 and J = ip 2 . The result is only a slight generalization of 



expressions that can be found in [66]. The generalization consists of allowing for a nonzero 
chemical potential (we also work in a different coordinate system.) These are then used to 
compute the on shell action and boundary terms that are needed in Section 4.2. 
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q 2 /j, 2 



2(2A_ - 1) 



r- 2 ) +/3r~ A+ ( 1 



2(2A 



+ 



a M „ r -3A- + 



(3A_ - 4) 



<j> = II- 


pr 1 + ^ 


1 = 1- 


a/3A_A+ 




6 


2 = 1- 


(m /2)r - ' 



q 2 fia 2 



hr 
fh 

where the two constants are: 
1 A A_(5A_-3) 



A_(2A_ - 1, 

a 2 A_ 



-2A_ i 

r + . . . 



4(2A_ - 1 



-2A_ i 

-r + 



a C A _ r -4A_ + 



(3A_ - 4) ' 



3A 



A 2 _ (26 A_ - 15) 



2A_ 



4(2A_ 



1} 



$(4A_-1) 32(2A_ - 1)(4A_ - 1) 



(D.i; 



(D.2) 



Following [6] the on shell Euclidean action can be shown to be a total derivative. The 
boundary terms give: 

S E = 0V(-2hh , f\ r=o °) (D.3) 

where V is the volume of space and /3 is the inverse temperature. The horizon does not give 
any contribution since / vanishes on the horizon. 

Following [6] the boundary terms needed in Section 4 are (they are all located at the 
boundary of AdS^j: 



S p ct = J d 3 x v/=^ (-2K + A^F^ - 2iPn»d^ + A B + m B ^ 2 + A B </> 4 ) 
= f3V h 2 (-f - Ah'h' l f + <f/<p - + Vfi^B + m B i? + \ B ^) 



lM''=oo (D _ 4) 

o 

(D.5) 

where is the metric on a fixed r slice as r — > oo. K is the trace of the extrinsic curvature 
and n M is the outward pointing unit normal to the boundary. By demanding that the Se + S^ 



is finite one can show that the constants in (D.4) are, 
Ab = 4, 



m B 



-A. 



A 



B 



3(3A 2 +4A) 
" 8(4A_ - 3) 



(D.6) 



Note that the ip A term in D.4 is only required when A_ < 3/4 exactly when C? 4 is relevant in 



alternative quantization. Plugging (D.I) into (D.3) and ( D.4[ ) one finds that the free energy 



is given by (4.6). 
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